Charm and bottom baryons in nonperturbative quark dynamics 
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We use the field correlator method in QCD to calculate the masses of Sc, He and 
recently observed T,fj, Ef, baryons and their orbital excitations. 
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^ ■ I. INTRODUCTION 

^ ! A comprehensive knowledge about the mass spectrum and spin splittings of heavy baryons 

is important for our understanding of quantum chromodynamics. The spectroscopy of c and 
P_|. b baryons has undergone a great renaissance in recent years. New results have been appearing 
in abundance as a result of improved experimental techniques including information on states 
r^, ' made of both light {u, d, s) and heavy (c, b) quarks [l]. Before 2007, the only one baryon 
with a b quark, the isospin-zero A^, was known. Now, we have the isospin one S;,jEj baryons 
O ■ and Sfo. The CDF Collaboration has seen the states and [2], while D0 ^ and CDF 



o 

00 

o 



J] have observed the S^. The masses of these states are summarized in Table [H 

On theoretical side there are many results on heavy baryon masses from different ap- 
proaches including a number of quark model variations , HQET [8:], sum rules [91] 
and lattice calculations Recently there have been several theoretical pa per s on the 
I masses of S;,, and Sf, based on the modelling the color hyper: ine interaction ll|, [l^, Il3 |. 
In the present paper we use the field correlator method (FCM) 1^ to calculate the masses 
of the S wave baryons containing c and b quarks and orbitally excited states that will be 
experimentally accessible in the future. The same exercise has been applied previously to 



the ordinary and cascade hyperons [15[. Using FCM, substantial work has been done in 



the heavy-light meson sector 16|] . However, so far only very few results have been reported 



for baryons. A further study of charmed and bottom baryons using FCM therefore seems 
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worthwhile. 

The dynamics of the ud pair plays a relevant role, being mainly responsible for the spin 
splitting in the strange sector. A similar contribution is expected for charmed and bottom 
baryons. Estimates of the one-pion exchange contribution to the baryon mass give -180 MeV 
both for A and A;,. Because our approach misses the chiral physics effects we calculate in 
this work the masses of the E and 5 states which are affected by the chiral dynamics only 
slightly. Note that the 6-baryons are structurally identical to the c-baryons: only a charmed 
quark is replaced by a beauty quark. Consequently, the analysis of b states and the results 
are only a variation of what is found for the charmed systems. 

II. THE FORMALISM 

FCM provides a promising formulation of the nonperturbative QCD that gives additional 
support of the quark model assumptions. The key ingredient of the FCM is the use of the 
auxiliary fields (AF) initially introduced in order to get rid of the square roots appearing in 
the relativistic Hamiltonian. Using the AF formalism allows to write a simple local form of 
the Effective Hamiltonian (EH) for the three quark system 

where Hq is the kinetic energy operator, V is the sum of the string potential and a one gluon 
exchange potential, rrii arc the bare quark masses, and /^j are the constant AF which are 
eventually treated as variational parameters. Such an approach allows one a very transparent 
interpretation of AF: starting from bare quark masses one arrives at the dynamical masses 
/li which appear due to the interaction and can be treated as constituent masses of quarks. 
The string potential considered in this work is 

Vriri, r2, rg) = armin, (2) 

where a is the string tension and r^j„ is the minimal length corresponding to the Y-shaped 
string configuration 

The mass Mb of a baryon is given by 



M = Mo + AEhf, 



(3) 
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where AEhf is the spin correction, 

^0 = E + f ) + Eoii^.) + C, (4) 

EQ^jii) being the energy eigenvalue of the Shrodinger operator Hq + V , and /Xj are defined 
from the minimum condition 

dMo{mi,Hi) ^ ^ 

For the light quarks (mj ^ -\/o") /Uj ~ v^(l + C'(q^s)) while for the heavy quarks (mj 
^/o^) fii ~ mj. In Eq. (jl]) C is the quark self-energy correction which is created by the color 



17j. This 



magnetic moment of a quark propagating through the vacuum background field 
correction adds an overall negative constant to the hadron masses: 

n 

where ri{t) is the known function [17} and 1/Xg is the gluonic correlation length. We use 
Xg = 1 GeV. 

Taking the approach implemented in the spin-independent masses can be obtained 
from (jl]). We solve the non-relativistic Schrodinger equation with the confining and Coulomb 
interactions by the hyperspherical method to determine the constituent quark masses /ij and 
the zero-order baryon masses Mq. Then we estimate HF splittings from the perturbative 
color-magnetic interaction with account of the wave function corrections. 

^ = E tt" < + 17 ■ V.))) (7) 

The first term in ([7]) is the standard color-magnetic interaction in QCD 18| , while the second 
term, proportional to the string tension a, was first derived in Ref. 19|. 

We use the basis in which a heavy quark is singled out as quark 3 but in which the 
light quarks are still antisymmetrized. The calculation of the spin matrix elements in (^^| 
is straightforward for J = 3/2, as the expectation value of each cTj cr^- is 1. For J = 
1/2 cri(T2 = lyCr^cTi = cr2 0-3 = —2 for and while for cti cr2 = —3, cr^cri = 

(72 (T3 = 0. 

The contact interaction in ([7]) requires the calculation of the 6 function expectation 
values. These contact probabilities were calculated using 3-body wave functions obtained 



4 



by a hyperspherical method. E.g. for L = the square of the baryon wave function at zero 
relative two quarks separation is 



{S{rij))L = o = /i-f ^7o, (8) 
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where 



oo 

2/ 



70 = / (9) 



is universal for all quark pairs, uo{x) is the hyperradial function normalized as 

u l{x) dx = 1, (10) 



with 
and 



= ft (r. - R^^r- = ^ rl, + 4, + ^ rl, (11) 



Hij = — , M = ni+H2 + fi3- 

Note that wave function corrections which influence the hyperfine splitting between the 
different baryons tend to affect 70 by only a few per cent: 70 = 0.1207 for nrac, 0.1197 for 
nsc, 0.1161 for nnb, and 0.1153 for nsb (in units GeV^^^). 
The second term in ([7]) is expressed in terms of the integrals 



{vij ■ Ki{\rij)) = — I u (x) 

TlAg 



00/2 \ 







dx, ^ = (12) 



iKi{Cl ^in'e cos' Odd 
V J 

III. THE RESULTS 

We employ some typical values of the string tension a and the strong coupling constant 
as that have been used for the description of the ground state baryons: a = 0.15 GeV^ 
and as = 0.39. We neglect the mass difference between u and d quarks, writing n to stand 
for either u or d. We use the current light quark masses m„ = 7 MeV and the (slightly 
updated) strange quark mass = 185 MeV found previously from the fit to Ds spectra 



16l |. However, our predictions need an additional input for the bare quark masses nic and 
rrib. These were fixed from the masses of Sc and S^, respectively, rric = 1359 MeV and 
rrib = 4712 MeV. 
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The result of the calculation of the 5* wave states is given in Table |TT1 In this Table we 
also present the dynamical quark masses fig and fig for various baryons. The latters 
are computed solely in terms of the bare quark masses, a and ag and marginally depend 
on a baryon. We also display the results obtained without the HF corrections. The baryon 
masses are for the isospin averaged states. The result show good agreement between data 
and theoretical predictions. 

In the Sq (with Q standing for either c or b) the light quarks are approximately in a 
state with S = 0, while another heavier state Sq is expected in which the light quarks 
mainly have S = 1. Both have total J = 1/2. The effect of S — S' mixing due to the spin- 
spin interaction is negligible 12|]. There is also a state Sq expected with total J = 3/2. 
The hyperfine splitting between S! and is found to be 69 MeV that agrees with the 
experimental value (~ 70 MeV) [20(], while the predicted mass difference — E![, = 26 
MeV agrees with the finding of Ref. However, our perturbative calculations do not 

reproduce the observed S'^ - Sc mass difference. The large hyperfine splitting between axial 
and scalar ns diquarks is usually described by the smeared 5-function that requires additional 
model-dependent assumptions about the structure of interquark forces. 

A similar calculations were performed for the P-wave orbitally-excited states, see Ta- 
ble IIVI Our basis states diagonalize the confinement problem with eigenfunctions that 
correspond to separate excitations of the light and heavy quarks (p- and A- excitations, 
respectively). Excitation of the A variable unlike excitation in p involves the excitation of 
the "odd" heavy quark. For states with one unit of orbital angular momentum between 
Q quark and the two light quarks we obtain M(Sc) = 2832 MeV, M(Sc) = 2867 MeV, 
M(Sfe) = 6132 MeV, and M(Sfe) = 6164 MeV, while the states with one units of orbital 
momentum between the two light quarks are typically ~ 100 MeV heavier. Note that zero 
order results of Table [IV] do not include the spin corrections (which are smaller than those 
for the S'-wave states) and the (negative) string corrections contributing into the masses of 
the orbitally excited baryons 2l|]. Our preliminary analysis of the latters shows that the 
string corrections tend to decrease the masses of the P-wave states by ~ 30 MeV. A more 
complete analysis will be given elsewhere. 
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IV. CONCLUSIONS 

We have calculated the masses of heavy baryons systematically using the FCM and 
the perturbative color-magnetic interaction. There are two main points in which we differ 
from other approaches to the same problem based on various relativistic Hamiltonians and 
equations with local potentials. The first point is that we do not introduce the constituent 
mass by hand. On the contrary, starting from the bare quark mass we arrive to the 
dynamical quark mass that appears due to the interaction. The second point is that for the 
first time we calculate the hyperfine splitting with account of the nonperturbative spin-spin 
forces between quarks in a baryon. We find our numerical results to be in agreement with 
experimental data and calculation in other approaches. 
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TABLE I: The masses of bottom baryons observed by CDF and D0 collaborations. 





Mass (MeV) 


Collaboration 




0oUo_2 gi^siai.j ± i. ( i^sysi. j 






5816;|:[!(stat.) ± 1.7(syst.) 




v*+ 


5829111 (stat.) ± 1.7(syst.) 


CDF[2j 




5837li:9(stat.) ± 1.7(syst.) 






5774± ll(stat.) ± 15(syst.) 


D0 [3] 




5793±2.5(stat.) ± 1.7(syst.) 


CDF [4] 



TABLE IL Heavy Baryons with L = 0. The values of the bare quark masses used in this calculation 
are nin = 7, = 185, nic = 1359, and m^, = 4712 MeV. The underlined masses have been 
used to fix m-c and mi,. The dynamical quark masses /ij are defined by Eq. ([5]). 



Bar yon 








Mo 






M 




470 




1455 


2479 


-19 


-6 


2454 




470 




1455 


2479 


30 


13 


2522 




476 


522 


1458 


2519 


-39 


-20 


2460 


E* 


509 
509 

514 


615 


4749 
4749 
4751 


5806 
5806 
5844 



+19 
-36 


2 

8 

-17 


5808 
5833 
5791 
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TABLE III: Masses of the heavy baryons from the present work and other approaches and the 
comparison with experimental data (in MeV). 







m 


m 


m 


m 


[9] 


this work 


exp 


S* 


2440 
2495 


2453 
2520 
2468 




2452 
2538 
2473 


2439 
2518 
2481 


2411+^5 


2454 
2522 
2460 


2454 ± 0.18 
2518.4 = 0.6 
2467.9 ± 0.4 


E* 

^6 


5795 
5805 


5820 
5850 
5810 


5824.2 
5840.0 
5805.7 


5847 
5871 
5788 


5805 
5834 
5812 


5809j:f^ 
5835+^y 
5780+61 


5808 
5833 
5791 


5808 
5829 
5774 ± 20 
5793 ± 3 
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TABLE IV: Heavy Baryons. L = 1. The A excitations involve the excitation of the "odd" quark 
(c for Ac, He or b for S^, H^), while the p excitations involve the excitation of the ud diquarks. The 
bare quark masses are the same as in Table HTl 



Baryon Lq, ^„ //^ Eq M 



nnc Ip 536 1452 1397 2920 

nnc 1a 495 1491 1377 2832 

nsc Ip 542 582 1455 1372 2954 

nsc 1a 497 544 1494 1353 2867 



nnb Ip 570 4746 1294 6240 

nnb 1a 540 4764 1234 6132 

nsb Ip 574 615 4748 1271 6272 

nsb 1a 542 588 4765 1211 6164 



